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Abstract. - We report the first experimental observation of reversals of a dynamo field generated 
in a laboratory experiment based on a turbulent flow of liquid sodium. The magnetic field ran- 
domly switches between two symmetric solutions B and — B. We observe a hierarchy of time scales 
similar to the Earth's magnetic field: the duration of the steady phases is widely distributed, but 
is always much longer than the time needed to switch polarity. In addition to reversals we report 
excursions. Both coincide with minima of the mechanical power driving the flow. Small changes 
in the flow driving parameters also reveal a large variety of dynamo regimes. 
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^ Dynamo action is the instability mechanism by which 
(~ — mechanical energy is partially converted into magnetic en- 
ergy by the motion of an electrically conducting fluid [1] . 
""^5 It is believed to be at the origin of the magnetic fields 
Oof planets and most astrophysical objects. One of the 
c/3 most striking features of the Earth's dynamo, revealed by 
£^>aleomagnetic studies [2], is the observation of irregular 
'"^reversals of the polarity of its dipolc field. This behaviour 
. . is allowed from the constitutive equations of magnetohy- 
. — drodynamics [1] and has been observed in numerical mod- 
^^ejs [3]. On the other hand, industrial dynamos routinely 
!_< generate currents and magnetic fields from mechanical mo- 
^tions. In these devices, pioneered by Siemens [4], the path 
of the electrical currents and the geometry of the (solid) 
rotors are completely prescribed. As it cannot be the 
case for planets and stars, experiments aimed at studying 
dynamos in the laboratory have evolved towards relax- 
ing these constraints. Solid rotor experiments [5] showed 
that a dynamo state could be reached with prescribed mo- 
tions but currents free to self-organize. A landmark was 
reached in 2000, when the experiments in Riga [6] and 
Karlsruhe [7] showed that fluid dynamos could be gener- 
ated by organizing favourable sodium flows, the electrical 
currents being again free to self-organize. For these exper- 
iments, the self-sustained dynamo fields had simple time 
dynamics (a steady field in Karlsruhe and an oscillatory 



field in Riga). No further dynamical evolution was ob- 
served. The search for more complex dynamics, such as 
exhibited by natural objects, has motivated most teams 
working on the dynamo problem to design experiments 
with less constrained flows and a higher level of turbu- 
lence [8]. The von Karman sodium experiment (VKS) is 
one of them. It has recently shown regimes where a sta- 
tistically stationary dynamo self-generates [9]. We report 
here the existence of other dynamical regimes and describe 
below the occurence of irregular reversals and excursions. 

Magnetic induction properties of turbulent swirling 
flows (VKS and others) have been widely studied ex- 
perimentally [8,10]. Dynamo action by these flows has 
also received strong support from numerical simula- 
tions [11, 12]. The VKS experimental set-up, shown 
in Figure 1, has been described in [9]. The fluid is 
liquid sodium, chosen for its high electrical conductivity 
(a = 10 7 ohm _1 m _1 ), low density (p = 930 kgm~ 3 ) and 
low melting point (98° C). The net volume of sodium is 
roughly 160 L. A turbulent von Karman flow is generated 
by two counter-rotating iron impellers (rotation frequen- 
cies Fi and F2). Its mean structure has the following 
characteristics: the fluid is ejected radially from the 
disks by centrifugal force and loops back towards the 
axis in the mid-plane between the impellers. A strong 
differential rotation is superimposed on this poloidal 



Berhanu et al. 




Fig. 1: The VKS2 set-up is designed to generate a dynamo 
flow in an electrically conducting fluid. The overall vessel is 
a copper cylinder of radius 289 mm and length 604 mm. The 
flow itself is confined within an inner copper cylinder (radius 
R — 206 mm, length 524 mm, thickness 5 mm), with sodium at 
rest between the inner and outer cylinders. An annulus of inner 
radius 175 mm (thickness 5 mm) is fixed along the inner cylin- 
der in the mid-plane between the disks. The counter-rotating 
iron impellers have radius 154.5 mm and are set 371 mm apart 
in the inner vessel; they are fitted with 8 curved blades of height 
h = 41.2 mm. Their rotation frequencies are independently ad- 
justable, up to 26 Hz. Magnetic measurements are made using 
a temperature controlled, 3D Hall probe mounted flush on the 
flow boundary, at the inner cylinder. 



flow, which generates a high shear in the mid-plane. The 
flow maximum driving power is 300 kW, and cooling 
is performed using an oil flow inside the copper walls 
of the vessel. It allows experimental runs at constant 
temperatures between 110°C and 160°C. The integral 
Reynolds numbers are defined as Re; = 2'KKR 2 Fi/v and 
take values up to 5 10 6 where v is the fluid viscosity and 
K = 0.6 is a coefficient that measures the efficiency of the 
driving impellers [12]. Corresponding magnetic Reynolds 
numbers, R m j = 2irKnoaR 2 Fi, up to 49 at 120°C are 
reached - /zo is the magnetic permeability of vacuum. 
The magnetic field is measured with local Hall probes 
inserted inside the fluid. 

When the impellers are operated at equal and opposite 
rotation rates F, a fully turbulent dynamo is observed 
when F is larger than about 17 Hz (R m = 31) [9]. The 
self-sustained magnetic field is statistically stationary with 
cither polarity in this case. In the experiment, the rota- 
tion rates (Fi,^) of the driving impellers can be inde- 
pendently adjusted and this gives an additional degree of 
freedom. Starting from a symmetric flow forcing, F\ = F2, 
one can progressively change the rotation frequency of one 
disk and explore regimes in which the faster disk imposes 
some kind of global rotation to the flow, a feature common 
to most natural dynamos. 

We show in Figure 2 a preliminary inspection of the 
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Fig. 2: (a) Preliminary inspection of dynamo regimes ob- 
served as the impeller rotation frequencies are independently 
set. Symbols: (o): statistically stationary dynamos, (+) no 
dynamo, i.e., magnetic field less than 10 gauss at the mea- 
surement location - and for measurement times longer than 
180 s. (★): dynamo with reversals, (b) Examples of the time 
variations of the main magnetic field component for rotation 
frequencies of the disks, 22 — 22 Hz (R m i = R m 2 = 42.5) and 
14 - 22 Hz (R m i = 27.5, R m2 = 43.5). Note that on measure- 
ment time scales of the order of 180 s, the regimes can depend 
on the path followed to reach them. 



parameter space accessible when the flow is driven with 
disks rotating at different speeds. As said above, only sta- 
tistically stationary dynamos are observed in the counter- 
rotating case (Figure 2b) . Another statistically stationary 
dynamo mode is observed when the frequency of one im- 
peller is increased from zero (say Fi), the other being kept 
fixed at 22 Hz, thus R m 2 in the range 42 — 43 depending on 
the sodium temperature (Figure 2c). Note however that 
its relative fluctuations are much smaller (compare Fig- 
ures 2b, c), an effect possibly ascribed to global rotation. 
This regime undergoes secondary bifurcations when the 
slower impeller frequency is increased further. In a small 
parameter range, AFi/Fi w 20%, a variety of dynamical 
regimes, oscillations, intermittent bursts (not shown), as 
well as dynamos with random reversals (Figure 3) are ob- 
served. We also find pockets of parameters for which we 
could not record the growth of a dynamo during 3-minutc 
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Fig. 3: Magnetic field measured inside the flow vessel, by a 3-dimensional Hall probe. No external magnetic field is applied, 
other than the ambient field, whose amplitude is about 0.2 gauss across the measurement volume. The temperature of the outer 
copper cylinder is T — 123 °C. (Main): Time evolution of all three magnetic field components. The main component (red) is 
the azimuthal one. Note that all components decay to zero at a reversal. The bottom graph shows synchronous recordings of 
the power driving the flow. (Right): detail of the time series of the main magnetic field and simultaneous power consumption 
(arrows mark the synchronous events). (Top): Chronos of the magnetic field orientation, white for a positive direction, black for 
the negative direction, for 2 successive recordings 900 and 1800 seconds long (separated by the shaded area, the first sequence 
corresponds to the main graph). In this regime, the von Karman flow is driven with counter-rotating disks at frequencies 
Fi = 16 Hz and F 2 = 22 Hz. 



long runs -corresponding to over 3000 forcing time scales. 

We now describe reversals of the magnetic field. In Fig- 
ure 3, we show a time series that corresponds to F\ = 
16 Hz and F 2 = 22 Hz. In this regime, the magnetic field 
reverses at irregular time intervals. All three components 
of the dynamo field switch polarity in perfect synchrony, 
so that B changes to — B. For each polarity, the amplitude 
of the magnetic field has strong fluctuations, with an rms 
fluctuation level of the order of 20% of the mean. This 
level of fluctuations is due to the very intense turbulence 
of the flow, as the kinetic Reynolds number exceeds 10 6 . 
Reversals occur randomly and have been followed for up 
to 45 minutes, i.e. 54000 characteristic time scales of the 
flow forcing. 

In the regime reported in Figure 3, the polarities do not 
have the same probability of observation. Phases with 
a positive polarity for the largest magnetic field compo- 
nent have on average longer duration ((T+) = 120 s) than 
phases with the opposite polarity ((T_) = 50 s). This 
asymmetry can be due to the ambient magnetic field. Note 



however that the amplitude of the magnetic field, that is 
much larger than the Earth field, is the same for both po- 
larities. Standard deviations are of the same order of mag- 
nitude as the mean values, although better statistics may 
be needed to fully converge these estimates. The mean 
duration of each reversal, r ~ 5 s, is longer than magne- 
tohydrodynamics time scales: the flow integral time scale 
is of the order of the inverse of the rotation frequencies, 
i.e. 0.05 s, and the ohmic diffusive time scale is roughly 
Trq ~ 0.4 s. Concerning the dynamics of field reversals, a 
natural question is related to the connection between B 
and — B in time. The equations of magnetohydrodynam- 
ics are symmetric under the transformation B to — B so 
that the selection of a polarity is a broken symmetry at the 
dynamo bifurcation threshold. The sequences of opposite 
polarities act as magnetic domains along the time axis, 
with Ising-type walls in-between them: the magnetic field 
vanishes during the polarity change rather than rotating 
as in a Bloch-type wall. 

One important discovery in these measurements is that 
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reversals of magnetic field are correlated with the global 
energy budget of the flow. The total power P{t) deliv- 
ered by the motors driving the flow fluctuates in time in 
a strongly asymmetric manner: the record shows short 
periods when P is much smaller than its average. They 
always coincide with large variations in the magnetic field, 
as shown in Figure 3. Either a reversal occurs, or the mag- 
netic field first decays and then grows again with its direc- 
tion unchanged. Similar sequences, called excursions [2], 
are observed in recordings of the Earth's magnetic field. 
The variation of power consumption during the weaken- 
ing of the magnetic field is in agreement with the power 
required to sustain a steady dynamo in the VKS2 experi- 
ment [9] (drops by over 20%, that is 20 kW out of 90 kW). 
However, we note that in other regions of the parameter 
space, different regimes also involve changes in polarity 
without noticeable modification of power. 

We have also observed that the trajectories connecting 
the symmetric states B and — B arc quite robust despite 
the strong turbulent fluctuations of the flow. This is 
displayed in Figure 4: the time evolution of reversals 
from down to up states can be neatly superimposed by 
shifting the origin of time such that B(t = 0) = for 
each reversal. Despite the asymmetry due to the Earth 
magnetic field, up-down reversals can be superimposed in 
a similar way on down-up ones if — B is plotted instead 
of B. For each reversal the amplitude of the field first 
decays exponentially. A decay rate of roughly 0.8 s _1 is 
obtained with a log-lin plot (not shown). After changing 
polarity, the field amplitude increases linearly and then 
displays an overshoot before reaching its statistically 
stationary state. 
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Fig. 4: Superimposition of 5 successive reversals from down 
to up polarity. For each of them the origin of time has been 
shifted such that it corresponds to B = 0. 

Further investigation of this regime will help address 
from an experimental perspective persistent questions 
about magnetic field reversals. Some of these concern 
the role of hydrodynamics and electromagnetic boundary 



conditions - both of them can be experimentally adjusted. 
Others are related to the dynamics of the magnetic re- 
versals. From inspection of paleomagnetic data, it has 
been proposed that reversing dynamos and non-reversing 
ones are metastable states in close proximity [13]. In 
gcodynamo simulations (convective dynamos in rapidly 
rotating spheres), the flow is often laminar and reversals 
have been associated to interaction between dipole and 
higher order modes, with the possibility of reversal 
precursor events [14]. Fields reversals have also been 
observed in turbulence driven numerical a 2 and aui dy- 
namos based on mean-field magnetohydrodynamics [15]. 
In these, the role of noise was found to be essential, 
together with the proximity of steady and oscillating 
states. In many cases, the existence of several dynamo 
regimes in a narrow region of parameter space has been 
considered as essential. Our experiment displays this 
feature: two different stationary dynamo modes bifurcate 
for F\ = i*2 and respectively F± ^ F%. Their interaction 
gives rise to a variety of different dynamical regimes in 
parameter space. This is a general feature for bifurcations 
of multiple codimension. The most striking aspect of our 
observation is that the low dimensional dynamics that 
result from the interaction of a few modes of the magnetic 
field is preserved despite strong fluctuations of the flow 
that generates the field. On average, the largest scales of 
the flow change more than a thousand times during each 
phase of given polarity. The large scale magnetic field 
cannot follow turbulent fluctuations and display features 
characteristic of low dimensional dynamical systems. 
Flows generating the magnetic fields of planets or stars 
involve dimensionless parameters (Reynolds and Eckman 
numbers) orders of magnitudes different from the ones 
of the present experiment and even further in the case 
of numerical simulations. However, a weak coupling 
between the large scale dynamics of their magnetic 
field and hydrodynamic fluctuations may explain why 
similar features are observed in some natural dynamos, 
in numerical modeling and in this experiment. 
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